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Abstract. We give a theoretical treatment of the interaction of electronic 
excitations (excitons) in biomolecules and quantum dots with the surrounding 
polar solvent. Significant quantum decoherence occurs due to the interaction 
of the electric dipolc moment of the solute with the fluctuating electric dipole 
moments of the individual molecules in the solvent. We introduce spin boson 
models which could be used to describe the effects of decoherence on the quantum 
dynamics of biomolecules which undergo light-induced conformational change and 
on biomolecules or quantum dots which are coupled by Forster resonant energy 
transfer. 
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1. Introduction 

Biomolecular systems contain a diverse range of optically active molecules that are 
crucial to biological function Important examples include retinal in rhodopsin 
which plays a role in vision, porphyrins and chlorophyls which play roles in 
photosynthesis photoactive yellow protein, blue copper proteins, and the green 
fluorescent protein (GFP) j3| which has become a powerful tool in molecular and cell 
biology. The different chromophores couple very differently to their environment which 
consists of the surrounding proteins and the solvent (water or aqueous electrolytes) 
leading to very different quantum dynamics 0]. Related questions arise concerning 
the spatial extent of the quantum coherence of excitations in spatially extended 
structures such as a-helices and /3-sheets in proteins and rings of porphyrin molecules 
in light harvesting complexes [H]. Thus, it is important to understand the coupling 
of electronic excitations in chromophores to the environment. This is not just of 
interest for understanding biological function, but also because these biomolecular 
systems could potentially be model systems to understand quantum decoherence at 
the nanoscale. This is because they have a complexity that cannot be fabricated 
even with the most advanced nanotechnology, but they can be modified and tuned 
using sophisticated biochemical techniques [HI- We also note that there are potential 
applications in biosensors with biomolecules coupling to quantum dots via resonant 
energy transfer 0, and with biomolecules acting as optically controllable elements in 
molecular electronics or as qubits in a quantum computer jlU) . 

Understanding the quantum dynamics of a system that is strongly coupled 
to its environment is a highly non-trivial, and sometimes controversial, problem. 
Furthermore, given the complexity and diversity of biomolecular systems it is desirable 
to find tractable model Hamiltonians which capture the essential physics and are 
amenable to rigorous analysis. The past few decades have seen monumental advances 
in combined quantum mechanical and molecular dynamics simulations, allowing the 
study of biomolecular systems of realistic size. See for example, the recent study 
of the quantum dynamics of retinal in the rhodopsin environment jllj . However, 
in spite of this progress it is unrealistic that in the next decade it will be possible 
to give a full treatment of the quantum dynamics including the quantum many- 
body effects associated with the strong interaction between the chromophore and 
its environment. Furthermore, simple models amenable analytical treatment can give 
insight into the essential physics involved with qualitative differences, especially when 
model parameters are extracted from quantum chemistry and molecular dynamics." 

We first consider the problem of the interaction of an electronic excitation 
(exciton) in a biomolecule (or a quantum dot) with the surrounding polar solvent. We 
look for the simplest possible model which can still describe the relevant physics, which 
allows us to determine the important parameters and perform analytical calculations. 
We show that this system can be described by an independent boson model where the 
environment is modelled as an infinite collection of harmonic oscillators with ohmic 
response. In the model we derive, the strength of the coupling to the environment is 
determined by the frequency dependent dielectric function of the polar solvent and 
by the difference between the electric dipole moments of the molecule in the ground 
and excited states. The relaxation rate of the polar molecules within the solvent 
determines the cutoff frequency for the ohmic response. 

We then introduce minimal models which can describe the effect of the solvent on 
the quantum dynamics associated with conformational changes at conical intersections 
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| 12l 1131 114| , and on excitations which are coupled via a mechanism such as resonance 
energy transfer. The latter is important in photosynthesis [2], is the basis of 
Fluoresence Resonance Energy Transfer (FRET) spectroscopy which is used to 
determine distances between chromophores in biomolecules |15l , and is the basis 
of new biosensors |7]. We have shown that each may be modelled by spin-boson 
models f3 El which exhibit rich many-body physics, and have been used to model 
systems such as the coupling between electron transfer and protein motion or a solvent 
^lElEOll^- The general spin-boson Hamiltonian is 

H = ieo-^ - ^AcTj; + ^ cjaO.'laa + ^cr^ ^ Ca{aa + al), (1) 

a a 

where ct^, (Tx are Pauli matrices describing the TLS, e is the energy separation of the 
two levels, A is the tunneling matrix element between the two levels, and the bath 
is modelled by boson operators acflj^ with frequencies uja and couplings Ca to the 
TLS. The independent boson model is with A = 0. This means that, unlike the 
spin-boson model, there are no transitions between the two quantum states. Only one 
environmental parameter, however, contributes to the system dynamics, the spectral 
density [TH1IT71 

JH = 47r^C2,5(c^-c^„). (2) 

a 

Leggett et al. jJSI defined a dimensionless coupling constant a which determines 
the extent of the quantum coherence. We show that for realistic biomolecules solvated 
in water this parameter will typically have values of order unity, which would prevent 
coherent Bloch oscillations for degenerate electronic levels. Such large values are in 
stark contrast to the Josephson junction qubits for which the effect of the Johnson 
noise in the electronic circuit is described by a spin boson model with values of a that 
are typically many orders of magnitude smaller |22| . 



2. Independent boson model for the interaction of a single chromophore 
with a solvent. 

We now show how an independent boson model can describe an electronic excitation in 
a biomolecule or quantum dot coupled to a solvent bath. We model the chromophore as 
a TLS with only a single active transition, of energy e. We assume the solute molecule 
has a permanent dipole moment in both the ground and excited state, given by fig 
and fie respectively. Describing the electronic ground and excited states by fermion 
creation operators and cj respectively, the solute part of the total Hamiltonian 
operator is 

Htls = ^e(4ce - clcg). (3) 

In the simplest picture, we can describe the solvation process by the well-known 
Onsager model ,23. .24) . This is a type of continuum model where the solvent is treated 
as a homogeneous dielectric, and has been used for both analytic and computational 
studies of solvation |25[I26) . The solute is treated as a point dipole which is surrounded 
by a spherical cage of polar solvent molecules with Onsager radius a (Figure^, which 
is typically the size of the solute molecule [331 IMl The cavity is assumed to 

be a vacuum, i.e., it has a dielectric constant = 1. The hydration shell of water 
molecules who have broken hydrogen bonds with their neighbours, and are now bonded 
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Figure 1. Onsager model of solvation. The chromophore is treated as a point 
dipole ft inside an empty cavity of radius a, which is surrounded by a polar solvent. 
The dipole polarises the solvent which in turn creates a an electric "reaction" field 
which acts back on the dipole, stabilising the solvated system. 



to the biomolecule, will serve to increase the effective cavity radius a. Further details 
(e.g., surrounding proteins) can be treated by adding shells with different dielectric 
constants around this cavity; this will be considered in a later work. 

The central dipole polarises the cage, which in turn produces an electric field 
inside the cavity, called the reaction field R{t). For the case of a uniform spherical 
cavity the reaction field is constant everywhere inside |2H1- This field acts back on 
the dipole, with interaction energy E = —t^{t) ■ R{t), typically lowering the total 
energy and hence forming a stable structure. Although the Onsager model does not 
include the microscopic details of the system, it does capture the essential physics of 
the solvation process jSH] and provided due care is taken to be aware of important 
microscopic effects (such as charge transfer processes) it should still provide valuable 
physical insight. 

Prior to the electronic transition, we assume the system to be in thermal 
equilibrium, and the reaction field will be parallel to the dipole fig. For simplicity, 
we assume that the direction of the dipole moment of the solute molecule in the 
ground and excited states points in the same direction. Defining = — f-g, 
Uz = c^Ce — cjcg, and identifying c]^Ce + cjcg with the 2x2 identity, 

Hsolute = ^eo-^, H.,nt = ^fT^(A^)i? + ^ (/ie + Hg)R. 

For a quantum mechanical model, we must quantize the reaction field. In the 
Heisenberg picture, we express the reaction field R{t) in its Fourier transform modes, 

R{t) = ^ e„ [a„e-''^=* + ale^^^'] , (4) 

a 

and quantise the coefficients so that and aj^ obey the boson commutation relations 
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We are therefore modelling the environment (solvent) as a bath of independent 
harmonic oscillators. The energy stored in the solvent cage is now expressed as 
^solvent = X]q '^'^Jt^.Q, and the full Hamiltonian for a single biomolecule coupled to 
the environment can then be written as 

H = ieo-^ + ^ UJaol^Ua + 

a 

a a 

where we define couplings Ma = ^e^A^ and = \ea{lJ-e + /^g)- 

This is the independent boson model POllS- We have not included a term, 
so that there can be no transitions between the two energy levels of the chromophore 
and the chromophore energy e (cr^) is a constant. Two conditions are needed to justify 
this. First, we consider only times after the initial excitation of the system (e.g., by a 
laser pulse) but much shorter than the radiative lifetime of the chromophore (typically 
1-lOns (32)). This means we do not include coupling of the TLS to the quantum and 
thermal fluctuations of the electric field which can induce optical transitions. At 
energies of the order of eV these fluctuations are much smaller than the dielectic 
fluctuations associated with the solvent. Secondly, for many electronic excitations we 
have assumed that the energy gap e (typically of order of leV) is significantly larger 
than the tunneling matrix element, and so its effects can be ignored on the timescales 
of interest. We defer discussion of it until below when we consider e becoming small 
near a conformational transition. 

The independent boson model can be used to investigate decoherence of systems 
coupled to a bath. Consider an initially coherent excitation, i.e., a coherent 
superposition of the ground and excited states of the chromophore. In the reduced 
density matrix of the chromophore, the population of each level is given by the diagonal 
elements, while the coherence is given by the off-diagonal elements. Thus the initial 
reduced density matrix of the chromophore has non-zero off-diagonal terms. Reference 
j83| shows that while the populations remain constant (as expected) the coherence 
decays exponentially, as given by Equations (18) and (19) and Figure 7 of 33 . Hence, 
after long times the reduced density matrix has no diagonal terms and so describes 
the chromophore in a mixed state |34) . 

The above model therefore describes the decoherence of an initially coherent 
excitation of a chromophore due to its environment (water, proteins), as well as 
spectral features such as the shifting and broadening of the absorption peak, as will 
be described below. 

3. Derivation of the spectral density function 

To complete the model we must specify the environmental coupling, i.e., the spectral 
function J(cl') from Following an approach similar to Caldeira and Leggett 
who considered a Josephson junction coupled to an electronic circuit, we relate J{uj) to 
the zero-temperature fluctuations in the uncoupled environment j21| . Noting that with 
no solute-solvent interaction {R{t)) = 0, we now examine the reaction field fluctuation 
correlation function S{t) IH^ defined as 

S{t) = i {R{t)R{0)) 6{t) = I (0|e'^*i?e~'-f^*i?|0) 6{t), (6) 
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where |0) is the ground state of the solvent harmonic osciUators, and H = Hsoivent — 
^^(jJadaCLa- We shall see that it is the imaginary part of the Fourier transform of 
S{t), which we denote f (w), that is relevant. We can write this as a sum over energy 
eigenstates [5^ . 

lm£iLu)^Try 6iLo-E,,)\{0\R\n)\\ (7) 

Expanding R into its normal modes again, we see that all terms (0|i?|n) vanish except 
for when a single oscillator is singly occupied. These states have energy i?„ = oja, and 
(0|i?|a) = Cq. Therefore, we see Im £{llj) ~ T^J^a ^a^i'^ ~ ^a), and using ((2Jl, we find 
Jiio) = {Afi)Hyji£{uj). (8) 

To calculate the reaction field fluctuations, we note that in the Onsager model |23[l37j . 
R{t) and the central dipole fi{t) are related by a linear response function x{^~t') ISHIi 
such that in Fourier space R{uj) = xi'-^) f^i'-^) ■ The static limit, where the solvent cage 
adjusts instantaneously, so that R(t) = Xs/z(t), is well known [281 123j . In reality, the 
solvent cage will lag behind a changing dipole due to dielectric friction |39[ I40|.and 
X(w) is given j2Z] by 

1 (9) 
' Aneoa^ 2eH + 1 ' 

where e{uj) is the frequency dependent dielectric constant of the solvent. We use the 
Debye formula [HI 112 

eM = eoo + ~ , (10) 

1 — tUJTD 

where Coo and are the high and low (static) frequency limits respectively, and tjj 
is the Debye relaxation time, the bulk reorientational relaxation time of the solvent 
dipoles J3^.44.- For water, these are given by Cs = 78.3, Coo = 4.21 and td = 8.2ps 

m 

We then apply the quantum fluctuation-dissipation relation 0^ for the 
imaginary part of x(w) 

X"{u;)^~m{l~e-P-)£{UJ), (11) 

which reduces at zero temperature to £{co) = 2ix"{co)- Note that this limit is well 
defined, even for w = 0. Using 11U|) and ijHJl, the spectral density is then given by 

^ ' 27reoa3 (2es + l)(2eoo + l)c^V| + 1' ^ ^ 

where te = ^ff^rn- 

This form for J(u;), which has a specific microscopic basis, is an important new 
result and could be used as an input into more phenomenological models such as 
the quantum Brownian model ^47 . In specific situations the magnitude of the J{lu) 
introduced here should also be compared to the J{uj) due to slow protein motion, as 
introduced in and ^I. 

We have thus shown that solute-solvent interaction in the Onsager picture can be 
modeled by a independent boson model with spectral density ifT^ . We note that J(ci;) 
has an implicit high frequency cut-off at Wc = related to the finite relaxation time 
of the solvent dipoles. Below the cut-off, J{oj) is approximately linear: J(w) = rjuj, 
where 

(A^i)^ 6(£, - eoo) 
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Thus, our model falls into the class of models defined in Ref. which below some 
high frequency cutoff have Ohmic dissipation, rj should correspond to the classical 
frictional coefficient. Indeed, rj is identical to the Nee-Zwanzig prediction of dielectric 
friction ^39. . 40 48, . 



4. Derivation of the optical absorption lineshape 



We can also use this model to derive an expression for the absorption lineshape a{uj) 
of the solute molecule in the presence of the solvent. In the gas phase, i.e., without 
the solvent, we would expect a sharp peak in the absorption spectrum at uj — e. In 
the presence of the solvent, however, this peak is shifted and broadened 0^ |^ ISU] 
by the solute-solvent dipole-dipole interactions. 

We use the results obtained for the independent boson model |30|. The resulting 
spectrum is 



a{Lu) 



-Re 



dt exp I it 



A- 



(t)]}, 



(14) 



where ^ is a normalisation factor (which we are not concerned with), and A = 
4 M^Ma ^ ^j-^-j jg -j-j-^g factor responsible for broadening the lineshape ^J. We 
express it ^1 07| using the spectral density |(5J : 



g{t) = J d.Ld 



J{i0) 



coth 



^1(1 



■ cos {cot}) + i sin {cut) 



Since J{uj) « for to 3> Wc, we can restrict the integral to u; < ujc, and note that 
for water at room temperature (300 K, m ~ 8 ps 02]) f^^ < P^c 1. We 
also make the short-time approximation j47| and Taylor expand the trigonometric 
functions to give g[t) = ^ksTXat^ + ^Xot, where we have defined the reorganisation 



energy Aq = -I / ^^^dw. We note also that because Ma 



M, we also have 



A = 2(a|^) ^0- Performing the Gaussian integral, we obtain the absorption spectrum 



a{Lu) = 



A 



where we have renormalised A, and 
(Aa^)' 6(e 



exp ■ 



UJ - {e- (^g/A^)Ao) 



An 



47reoa3(2e« + l)(2eoo + l)' ^^^^ 

If we set Eoo = 1 then the resulting lineshape is precisely that of where the 
absorption spectrum was obtained classically by assuming a Boltzmann distribution 
of the environment. Thus, our quantum model gives the correct classical limit. 



5. Spin-boson models of conformational change through conical 
intersections 

This model of solvent interactions can be extended to other situations. A number of 
important photochemical reactions in biomolecules involve a conformational change 
that initiates a signal. Examples include the cis/ trans isomerisation of retinal |13) and 
the photoactive yellow protein (PYP) j51| . The former drives the proton pumping (in 
rhodopsin) against an external pH gradient that leads to the signal responsible for 




Figure 2. A plot of two energy surfaces intersecting in a conical intersection. 
The vertical axis is energy, and the horizontal axes are qi and 52 representing the 
position co-ordinates of two classical harmonic oscillators. Scales are in arbitrary 
units. 



vision. The latter causes a conformational transition of the surrounding PYP protein. 
Only recently has it been shown [I'M |22I that these isomerisations occur through a 
"conical intersection" |I2j, between the potential energy surfaces of the ground and 
excited states. 

A minimal model for this requires two nuclear co-ordinates (degrees of freedom) 
|52|. which form a two-dimensional simple harmonic oscillator, corresponding to the 
classical "cone shaped" potential surfaces [21 ■ The simplest model, with classical 
nuclear co-ordinates, which can describe a conical intersection is 

H ^ecr^ + Acr^ + -^{pl + mlujlql) + -^{pl + mlujlql) + cjaCzqa + .9f,crx9f,, (16) 
Zma Znib 

where ga,b and Pa,b represent the position and momentum respectively of the two 
classical oscillators, cja.b and ma.b are the frequencies and masses, respectively, of the 
oscillators, and ga^b is the coupling of each oscillator to the chromophore. This model 
includes a coupling term Aa^ between the two states, where A is the matrix element 
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responsible for the non-adiabatic mixing of the ground and excited states |12j. The 
spectrum is 

E± = ±\j{e + gaqaf + {^ + 9bqbf + \{ojIpI + ql+ uj^^pI + 9h) (17) 

which clearly has a conical intersection when the first term vanishes (Figure |2Jl- We 
note that for e = this corresponds to the E ® e Jahn- Teller model jHSj • Hahn and 
Stock |14| recently studied such a minimal model for retinal in rhodopsin but did not 
include the effect of the solvent, which can play a pronounced role on the topology of 
the conical intersection [^. The generalisation of this model to include the solvent 
and quantum motion of the nuclear coordinates is 

H =]^^(Jz + Actj, + ct^ ^ Ma (al + aa) + X! ^""q"" + (1^) 



WaO^a + LJbb^b +ga\ -^^<yz{a + a)) + \ -^^gbiJxih + b'^), (19) 

V maLda V mbUJb 

where a) and 6^ are creation operators for the vibrations for the two degrees of freedom 
associated with the conical intersection. 

An important question about the quantum dynamics of the chromophore retinal 
concerns the explanation for the following observation jS^]- The speed and quantum 
efficiency of the conformational change is much greater when the retinal is in the 
rhodopsin environment rather than in a solvent. This model could be used to address 
this question by performing the following (non-trivial) calculation. Suppose following 
optical excitation the TLS is in the excited state and the bath and vibrational modes 
are in the equilibrium state associated with the ground state of the TLS. We want to 
calculate the branching ratio between the two possible final states of the system: to 
the ground electronic state and the excited electronic state and the associated phonon 
equilibrium state. Transition to the former requires that the non-adiabatic coupling 
term induce a transition to the ground state by a mechanism similar to Landau-Zener 
tunneling. However, in the presence of strong decoherence this tunneling will be 
suppressed. 

We now consider the magnitude and effect of the coupling of the bath to the 
TLS near the conical intersection. In the case of Ohmic dissipation, we can define a 
dimensionless coupling a — rj/2nT%. For the spin-boson model ^ with A ^ 'hLOc^ if 
a > ^ and e = there will always be incoherent relaxation [IB) , while for a > 1 and 
small A the system is localised at its initial state. For e = and a < i, coherent 
oscillations between the two states will occur. These will be inhibited by increasing 
the bias e and the temperature TH but enhanced by values of A much larger than 

We can make an estimate for a in water, using the constants given earlier which 
completely describe the solvent. The two solute unknowns are the cavity size and 
change in dipole moment. If we measure these in angstroms and Debye respectively 
{ID — 3.3 X 10^^°C • to) then a = 22 ^^^^ . For a typical chromophore with radius 
between 5 — SA, a dipole moment change of just 1 — 5D is sufficient to make a > 1, 
and this condition seems likely to be met for most small molecules IS^l Ch. 6]. In 
contrast, we note that for GFP, quantum chemistry calculations suggest A/i is small 

The solvent is therefore a large potential source of decoherence which destroys 
any quantum coherent oscillations of the electronic excitation. 
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2a — • — 2b 2a — I — 2b 

Figure 3. Schematic diagram of the non-radiative transfer of an excitation from 
molecule a to molecule b. Initially, molecule a is in the excited (2a) state, while 
molecule b is in the ground state (16). When coupled by fluorescence resonant 
energy transfer, the first molecule is de-excited to the ground state la, while the 
second molecule transitions to the excited state 2b. No photons are exchanged 
during the process. 



6. Fluorescence resonant energy transfer 

Another important process in biomolecules is coherent transfer of excitons. An 
example is fluorescence resonance energy transfer (FRET) , as shown in Figure |31 
This occurs in photosynthesis T and is the basis of FRET spectroscopy [151 and 
a proposal for quantum computing based on biomolecules |10| . If 1 and 2 denote the 
two molecules, the coupling can be described by adding a term to the Hamiltonian 

Ho = ^ {a\al + ala\) = J12 + ay^) , (20) 

where = if^y- J12 describes the FRET coupling and scales like where R 

is the separation of the molecules |S7] . 

What about the effect of the solvent? We assume that the molecules are 
sufficiently far apart that the reaction fields acting on the two biomolecules are 
uncorrelated. This assumption can be tested by looking at the wavevector dependence 
of the dielectric constant e(w, g); this will be weak. We anticipate that the 
corresponding length scale is of the order of the separation of the polar molecules in the 
solvent which in turn is much larger than i?, the separation of the biomolecules. Thus, 
we model the reaction fields by two independent collections of harmonic oscillators. 
The relevant Hamiltonian is then 

H = ^eial + UJaaioa + ct] ^ ga{aa + a],) + 



a a 

Ji2(a>2+^V2), (21) 



and J{uj), J{uj) are as before, and ei, €2 are the energy gaps of the two biomolecules. 
In the Ohmic regime, this will have similarities to two Kondo impurities coupled by an 
XX interaction and in an external magnetic field [SS]. We have studied this problem 
in more detail in a separate article 59 . There, we show that in most cases the above 
Hamiltonian may be reduced to a single spin-boson model. This allows us to give 
quantitative criteria, in terms of experimentally measurable system parameters, that 
are necessary for coherent Bloch oscillations of excitons between the chromophores, 
and suggest how these could be observed experimentally through FRET spectroscopy. 

The above Hamiltonian (|21|l has a natural generalisation to many coupled 
biomolecules, such as the rings of chlorophyll molecules in a light harvesting complex 
in purple bacteria or in a segment of base pairs in double stranded DNA jHO]. It 
will allow calculation of the localisation (the extent of spatial coherence) of excitons 
in such extended systems. 
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7. Summary 

The major new result of this work is that we have shown that a chromophore and 
its environment may be modeUed by an independent boson model, and derived the 
appropriate spectral density, which can be expressed in terms of experimentally 
obtainable parameters. We have estimated the critical coupling parameter a, and 
found that quantum states are rapidly decohered by the interaction of the chromophore 
dipole with its surrounding polar environment. However, as can be seen from equation 
H13|l . the coupling to the environment is strongly dependent on the cavity radius, 
and further work must be done to determine the effects of surrounding proteins and 
solvation shells which might push the bulk solvent away from the chromophore and 
reduce the decoherence. Finally, we have introduced spin-boson models which could 
be used to model conformational change around a conical intersection and molecule 
coupled by fluorescence resonance energy transfer. These models, although simple, 
exhibit rich physics and may help to provide insight into the workings of these complex 
systems. 
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